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Abstract 


The report presents the application of fixed lag smoothing algorithms 
to the problem of estimation of the phase and frequency of a sinusoidal 
carrier received in the presence of process noise and additive observation 
noise. A suboptimal structure consists of a phase-locked loop (PLL) 
followed by a post-loop correction to the phase and frequency estimates. 
When the PLL is operating under high signal-to-noise ratio, the phase 
detector is approximately linear, and the smoother equations then 
correspond to the optimal linear equations for an equivalent linear signal 
model. The performance of such a smoother can be predicted by linear 
filtering theory. However, if the PLL is operating near the threshold 
region of the signal-to-noise ratio, the phase detector cannot be assumed 
to be linear. Then the actual performance of the smoother can 
significantly differ from that predicted by linear theory. In this report 
we present both the theoretical and simulated performance of such smoothers 
derived on the basis of various models for the phase and frequency processes. 
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I. Introduction 

The derivation of optimum receivers through modern estimation 
techniques has been proposed by various researchers [1-12] (see also their 
references). In references [4,5] optimum zero lag receivers have been 
derived on the basis of linear Kalman filtering theory [6] for linear 
measurement schemes. The nonlinear measurement situations which are of 
interest here have been studied in [7], wherein, on the basis of nonlinear 
filters of [8], suboptimal nonlinear zero lag receivers have been derived 
for the demodulation of angle modulated signals. In [9] the techniques of 
[7,8] have been extended to design suboptimum nonlinear fixed lag 
smoothers for phase estimation (see also [10-14]). The solution of the 
optimum nonlinear filtering and smoothing problem is, of course, 
intractable. 

Whether derived from linear or nonlinear theory, the smoother 
structure consists of a phase-locked loop (PLL) followed by a post-loop 
correction to the phase and frequency estimates. In this report we study 
the application of linear and nonlinear smoothers to the estimation of 
phase and frequency of the sinusoid. We show that for this case, the 
suboptimum nonlinear smoother equations derived from [7*8,9] are not 
substantially different from the optimum linear smoother equations. In 
addition, simulations show that the difference in performance of the 
nonlinear and linear smoothers is small. 

In this report we evaluate the performance of smoothers both by 
analysis and simulations. For high signal-to-noise ratio conditions when 
the phase detector can be assumed to be linear, the results predicted from 
linear theory are in close conformity with those obtained from simulations. 
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However, when the phase detector is highly nonlinear, the linear theory is 
inadequate to predict performance even for linear smoothers. 

The performance of the smoother with a nonlinear phase detector is 
evaluated by simulations. When compared to the case of a linear phase 
detector the performance of the smoother can be substantially different. 
The difference can be much more pronounced when the process noise is 
present, compared to the case when only the observation noise is present. 
The simulation examples indicate that in the absence of process noise, 
although there is a significant performance degradation due to nonlinearity 
(about 1 dB when operating in the loop SNR of about 5 dB), there is no 
threshold observed in the smoothing error covariance in this region. In 
contrast to this, there is a pronounced threshold in the smoother 
performance when the process noise is present and the inverse of filter 
phase error variance is below 7.5 dB. 

We also evaluate the smoother performance when the process noise is 
reduced in magnitude. That is, the smoother/ filter solutions are based on 
a relatively high process noise, but in the simulation the actual variance 
of the process noise used is lower or zero. This is of interest because, 
in many practical applications, the process noise statistics are not 
precisely known, and are therefore deliberately overestimated. 

In Section II we present the signal model and the suboptimum smoother 
equations as derived from the stochastic optimization theory. With a minor 
modification these equations are shown to be the linear Kalman filter 
equations for an appropriate linear signal model applicable when the phase 
detector can be assumed to be linear. This section also presents a simple 
finite impulse response (FIR) filter implementation of the smoother. Under 
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steady-state conditions, the coefficients of the FIR filter are constant 
and can be precomputed. 

Section III presents three specific subcases of the signal model that 
are commonly used in the present application including those involving the 
phase and frequency estimation for communication and navigation systems. 
In this section, assuming the phase detector to be linear, closed-form 
expressions are derived for the filter and smoother transfer functions and 
their performance indices. 

Section IV contains the simulation results for both linear and 
nonlinear phase detectors and compares the simulation results to those 
obtained from theory. 

Finally, some concluding remarks are made in the last section. 
Detailed derivations are contained in Appendices A to C. Appendix D also 
contains the flow chart for the computer program developed for the above 
simulations and computations. 
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II. Signal Model and Smoother Equations 

We consider the problem of estimating the phase process 0(k) from the 
sampled version of the received carrier signal y(k), i.e., 

y(k) = A VT Sin (w e t k + 0(k)) + v(k) (1) 

where t k is the kth sampling time, w c is the known carrier frequency and 
v(k) the observation noise is the sampled version of a narrow band zero 
mean white Gaussian noise process v(t). Furthermore, the phase process 

0(k) is modeled as 


e (k) = Bfi'x(k) ; e= CIO 0] 


( 2 ) 


x(k+1 ) = $ x(k) + w(k) 

In (2), 6 is the phase constant, x(k) is the state vector of dimension n, 
$ is an (n x n) matrix and w(k) is zero mean white Gaussian noise process 
independent of {v(k)}. Thus 

E[v(k)] = 0 , E[w(k)] = 0 

E[v 2 (k)] = R ; E[w(k)w T (k) ] = Q ; E[v(k) w(j)] = 0 


Three possible values of n are considered in the sequel corresponding to 
first, second and third order smoothers. The specific selection of the 
matrices Q and $ is considered later. 

Applying the imbedding approach of [14], in reference [93* the 
suboptimal nonlinear equations have been derived from the nonlinear filter 
equations of [8]. These equations are presented in appendix A. 
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As shown in appendix A f when 2<a Q t and higher order harmonic terms are 
ignored, the smoother equations reduce to the following 

x (kfl) = $£ 0 (k/k) + K 0 (k+l)n(k+l) 
x^(k+l/k+l) = x^ ^(k/k) + K^(k+l)n(k+l) 
n(k+l) = /Ty(k+l)Cos(0(k+l)) 

0(k+l) = U) + 3c'x(k+l/k) (3) 

c k+l 

x^(k) = x(k-i) ; i=0,...L 

In (3) x ± (k/J) denotes the estimate of x^k) on the basis of observations 
up to time j and the gain vectors and the cross covariance matrices 
P i0 (k/j) = E{x(k-i/j) x T (k/j)} are given by 

K^k+1) = A 3P 10 (k+1/k) 8 S-^k+l) ; 0 < i < L (4) 

P i0 (k+l/k+l) = P io (k+l/k) - P. 0 (k+l/k)(AB8)(ASiO' P 00 (k+l/k)S _1 (k+l) (5a) 


P. 0 (k+l/k) = p i_ 1>0 ( k / k ) 5 0 < i < L 


(5b) 


p a = S 2 p n i’ 1 (6) 


<!> 


00 


<P 

where R(k) = R(k)/A iJ and the smoother error covariance matrix P ii (k+1/k+1) 
is 

P i; .(k+l/k+l) = P..(k+l/k) - P. Q (k+l/k)(AB8) S _1 (k+1)(A6^)" P" 1Q (k+l/k) 
P..(k+l/k) = P^.^k/k) 


(7) 


A. Rapprochement with Linear Theory: 

Representing the bandpass additive noise v(k) in terms of its baseband 
quadrature components v^(k) and Vq(k) as, 
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v(k) - (k)Cos(o) c t^ + $£/x(k) ) + v (k)Sin(w t + B£"x(k)) 

i C K 

and ignoring the 2a Q term, the phase detector output is given by 

n(k+l) = A Sin(erx(k+l/k)) + V F==- ^v i (k+l)Cos(3£'5c(k+l/k)) 

- v q (k+l)Sin(B£'x(k+l/k))J 

x(k+l/k) = x(k+l) - x(k+l/k) 

For small estimation error x(k+1/k), 

r|(k+l) = BA£^x(k+l/k) + ^===. V i (k+1) ^ 

The n(k+1) given by (8) above is precisely the one-step ahead prediction 
error (innovation) for the following linear model 

y(k+1 ) = BA£"x(k+l) + — v.(k+l) (9) 

VT 1 

It is easily verified that equations (3-7) reduce to the linear 
optimal smoother equations for the model (2,9). 

B. Smoother Implementation: 

Figure 1 illustrates the smoother implementation. If the various 
gains are replaced by their respective steady state values, the smoother 
consists of a sampled data phase-locked loop followed by a post-loop 
correction to the filtered estimates. There is an equivalent and somewhat 
simpler implementation to that of Figure 1. To derive this, one notes from 
figure 1, or the equations (3)» that 

J(k+l-L/k+l) - i 0 (k+l-L/k+l-L) + K L (k+l)n(k+l) + 

+. . .+K 1 (k+2-L)n(k+2-L) (10) 
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Under steady state conditions it is easily seen that 


K. = (6A S~ 1 P F )y 1 1 , P F = lim P 0Q (k/k) 

k -*•“ 

Y-1-8AP S ; Pp = lim Poo(k+l/k) (•)■]) 

k-* 00 

Substitution of (11) in (10) results in the implementation shown in Figure 
2. The finite impulse response filter (FIR) has output £ (k+1 ) r elated to 
its input n (k+1 ) by 

e(k+l) = n(k+l) + y~ 1 Tl( k ' ) +•••+ Y ^ L_1 ^ ri( k +2-L) (12) 

and ip in Figure 2 denotes the vector 


^ = 6A S -1 P f Y L_1 £ 
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III. Derivation of Steady State Filter/Smoother Transfer Functions and 

Performance Expressions 

In this section we consider three specific cases of model (2). The 
resulting filter/smoother configurations are termed first order, second 
order and third order respectively. By replacing various gains and 
matrices by their steady-state values in equations (3-7), these difference 
equations are replaced by algebraic equations and may be solved explicitly 
for the steady state values of the filter error covariance matrix Pp, the 
prediction error covariance Pp, smoother error covariance Pg» etc., as in 
[15-17]. In the following, we consider the measurement model (9) 
(justified above) and replace g kl' by h’ for notational convenience. 

A. First-Order Filter/Smoother : 

In this case x is a scalar denoting the unknown phase, $ =1 and the 
scalar H may be assumed to be 1 without loss of any generality. The 
derivations for P p , P p , etc., are simple for this case, one may refer to 
[16] for example. Substitution of these in equations (3) also yields 
transfer functions of the steady state filters and smoothers as in [15]. 


B. Second-Order Fil ter/ Smoother : 
For this case, 


$ 


1 

0 


T 

1 


H • = BA [1 0] 


(13) 


The state vector consists of phase and frequency to be estimated. In [15- 
17], expressions for Pp and Pp are derived for a very specific Q matrix. 
Here, we also derive such expressions for any general matrix Q. This 
derivation is contained in Appendix B. Substitution of these and other 
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related parameters in (3) yields the various transfer functions associated 
with the steady state filters and smoothers. The detailed derivation is 
contained in Appendix C. The result for steady state smoother error 
covariance matrix is also given in Appendix B. 


C. Third-Order Filter/Smoother: 

For this case, 

' 1 T T 2 /2" 

_ • 

0 1 T 

0 0 1 


a' = $a [i 


0 0 ] 


(14) 


The state vector consists of phase, frequency and the frequency 
derivative to be estimated. The expressions for steady-state values for 
various covariance matrices may be derived as for the case of first order 
and second order smoothers/filters. Alternatively, thes'e are also 
obtainable from the simulation of the algorithm. Substitution of resulting 
steady-state values of Kalman gain in (3) results in the transfer 
function of the steady-state filter/smoother as shown in Appendix C. 

The above derivations are based on the assumption of a linear phase 
detector. The performance predicted on the basis of these expressions 
(Appendix B,C) is compared with simulations in the next section. As would 
be observed there, under the assumption of a linear phase detector, the 
simulation results are in close conformity with those predicted from 
theory. 
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IV. Simulation Results for Linear and Nonlinear Cases 


Appendix D contains the flow chart of the simulation package developed 
on the VAX system for evaluating the performance of the filters/smoothers 
by simulations. The program includes both time domain and frequency domain 
analysis. Thus from the steady-state filter covariance matrices obtained 
from simulations, the program can also provide coefficients of various 
transfer functions associated with filters/ smoothers for various delays and 
their two-sided noise bandwidths. In the following the simulation results 
are presented for the performance of smoothers. First we summarize the 
main points of these results in section V. This is then followed by a more 
detailed presentation in section VI. 

V. Simulation Results-Summary 

The simulation results obtained for a second order smoother can be 
broadly categorized into linear and nonlinear cases as follows. 

A. Linear Cases 

The first part of the simulations/computations evaluates the 
performance of the filter and smoother with a linear phase detector. Here 
we verify the results predicted from theory with those obtained from the 
real time covariance matrices and the simulations. The two sets of results 
agree within the bounds of statistical errors. Moreover from these 
performance plots, simple expressions for the filter/smoother noise 
bandwidth are also derived in terms of process noise covariance and the 
sampling period. This is done for both the phase and frequency estimation 
error variances. We consider three specific cases for the process noise 
Variances. The first case corresponds to zero process noise. The second 
corresponds to nonzero process noise with its covariance known and in the 



third case we consider the situation where the process noise variance is 
overestimated. The results show that the frequency estimation error 
variance is a much stronger function of the process noise than is the phase 
estimation error variance. Thus the noise bandwidth of the frequency 
estimator varies as some power £ of noise bandwidth of the phase estimator. 
For the simulation examples, £ varies between 3 to 6 for the sampling 
period T between 0.01s and Is. We also establish that the asymptotic 
smoother improvement over the filter performance is about 5.6 dB for both 
the phase and frequency estimation. Another interesting result of 
simulations is that for phase estimation alone, a normalized (with respect 
to the dominant filter time constant) smoother delay of 1 is adequate to 
achieve asymptotic performance for the simulation examples. To improve the 
frequency estimates as well, the normalized delay is approximately 2. 

B. Nonlinear Cases 

The second part of the simulations evaluates the filter and smoother 
performance when the phase detector nonlinearity is taken into 
consideration. Here we observe that when the actual simulated process 
noise variance is zero but the filter/smoother is designed on the basis of 
nonzero process noise variance, the smoother still provides an improvement 
of about 5.1 dB or more when the loop SNR (linear theory) exceeds 6 dB. 
However for this range of loop SNR, the filter itself may result in a 
degradation of 1.5 dB or less due to the nonlinearity. Thus overall with 
the smoother and for a given value of noise bandwidth (which may be 
constrained from other considerations) the receiver may operate with more 
than 3«5 dB reduction in the carrier power to noise spectral density ratio, 
if the required phase error variance is about 0.1 or smaller. A similar 


result is shown to be valid when the simulated process noise variance is 
nonzero. For the simulation example with loop SNR J> 7.5 dB, similar 
improvements are obtained. However, for loop SNR less than 7 dB, the 
smoother exhibits a very marked threshold effect in that the smoother 
improvement becomes small with the reduction in the loop SNR. 
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VI. Simulation Results-De tailed Presentation 


In the following the simulation results obtained for the second-order 
case are presented in some detail. To be concrete we use the following 
often used model for the Q matrix 


Q = 


T ^ /3 
T/2 


T/2 

1 




(15) 


One advantage of using the above Q is that the performance of the 
filter /smoother is then the function of only three parameters viz, o 

Si 

2 2 

a-, and T, where a- = a /2 and denotes the noise variance of 

v v v 

v^(k)/V2~ in the baseband model (9)« We present the smoother/f il ter 
performance in terms of both the phase estimation error and frequency 
estimation error variances. 


A. Phase Tracking Performance 

First we present the results from both analysis and simulations for 
the case when the phase detector nonlinearity is ignored. Both the optimal 
filter and smoother performances are analyzed in the following. 

1. Optimal Filter Performance 

Here we present the performance of the optimal filter (assuming linear 
signal model) so as to relate these parameters to the two-sided normalized 
loop noise bandwidth 2Bp, (subscript p denotes phase) a commonly used 
parameter in the design of phase-locked loops. In the simulations, A and B 
have both been normalized to 1. Thus in terms of the closed-loop transfer 
function matrix Gp(Z) of (C8) (with Kg set equal to zero), the parameter 
2B p is given by, 
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2Bp - -i- | S,p (Z)G F,p <rl) -7- (16) 

r 

where Gp >p represents the first component of the transfer function matrix 

Gp(z), and r is some appropriate contour of integration. Figure 3 plots Bp 

2 2 

as calculated from (16) as a function of (a /a- ) for three different 

a v 

values of the sampling period T viz. 0.01s, 0.1s and Is. From these graphs 
it is readily seen that 


2 2 - 11 

2B .67(a /a- ) T 

p a v 

, 9 .226 

2B 35 i.01(a /O- ) T 

p a v 

2 2 ' 

2B 35 1.05(0 /O- ) T 

p a v 


T = Is, (o 2 /o- 2 ) < 10 
a v 

T = 0.1s 
T = .01s 


(17) 


Equation (17) shows that the relation between the real two-sided bandwidth 

2 2 

2B./T and the parameter (o /o- ) (not a function of T) is also dependent 
P a v 

upon the sampling period T. The last relation in (7) may be taken to be 
the asymptotic relation for (2B p /T) as T — ► 0. Figure 3 also includes 
the normalized loop noise bandwidth of the smoother as calculated from (16) 
with Gp^p replaced by the first component of the transfer function matrix 
given by (C14-C16) with Kg=0. Denoted by 2Bp,g this normalized bandwidth 
is given by, 

2 2 .234 

2B P ,S ~ - 28(CT a /a v > T 5 T = 0.1s (18a) 


Comparison with the filter bandwidth of (17) yields 


(B/B „) — 3.6(ct 2 /o- 2 ) 
P P»S a v 


-.008 


, T = 0.1s 


(18b) 


Figure 4 plots the real two-sided noise bandwidth of the filter and 
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smoother denoted by 213^ and 213^ g respectively where B^ p = B p /T and B Lp ^g 

= Bp^g/T. The figure also plots the normalized value of the estimates of 

phase error variances both for the filter and smoother. These variances 

are the steady-state values of the (1,1) components of the actual filter 

error covariance Pp and the smoother error covariance matrix Pg 

2 

respectively and have been obtained by simulation, with a set equal to 

a J 

2 2 

zero. The term a t denotes the variance of the samples of the second 

ab 

component of the process noise w(k) actually used in simulations. The term 

o_ 2 in (15) thus represents an estimate of a 2 for the purposes of 

a a,S 

2 

filter /smoother equations (design) and thus in general o may be different 

a 

2 

from ° a s • In case the process noise covariance matrix is known 

2 2 

precisely, then of course a = a c . As is evident from Figure 4, 

a a j d 

2B Lp - P F c 1 , 1 ) / (a— 2 T) ; 2B Lp>s = P s (l,l)/(a- 2 T) (19) 


where ~ denotes estimate. 

In Figure 5 we plot the normalized phase error variance for both the 
filter and smoother as obtained from the recursive solutions of equations 
(4-6). From the figure approximate expressions for these terms may be 
written as, 


P (l,l)/(a- T) a 

t V 


2 2 
.75(a /a- ) 

a v 

1.32 (c 2 /a- 2 ) 
a v 


.08 


.22 


1.4 (a 2 /a- 2 ) 

a v 


.25 


T = Is , (a 2 /a- 2 ) < 10 
a v 

T = 0.1s 
T = .01s 


(20a) 
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.237 


(20b) 


P (1 , 1) /(a- 2 T ) 2 .365 (a V- 2 ) 

o v a v 


Comparing (20) and ( 17 » 1 8) one observes that provided an optimum 

2 

filter or smoother is used corresponding to the process noise variance o a , 

the maximum degradation of the phase error variance (for a specific loop 

2 2 

noise bandwidth) due to nonzero process noise variance c < a 

3 3 

is only about 1.34 (1.25 dB) compared to the case of o 2 = 0, and this is 

2 

almost independent of the variance cr . 

Figure 6 provides a direct comparison of the filter performance with 

2 

and without a _ . This is in view of the fact that under the assumption 
3 9 S 

2 

of linear phase detection, and with o c = 0, the phase estimation error 
variance P(1,1) s (NqBl/P c ), where P e is the received carrier power (P c s 
A and is equal to 1 for simulations) and N Q is the one-sided power 
spectral of the noise v(t) in (1). With a- 2 = a 2 /2 = N Q /(2T), we have 

Pp(1,1)/(o- 2 T) = 2Pp(1,1)/Ng = 28^ establishing the equivalence of the two 

2 

performance measures under the assumption of a ag = 0. Thus, in figure 6, 

2 2 

the plots of 2 Bl correspond to Pp(1,1)/(cf- T) with a a s = °» while those of 

2 2 2 

Pp(1,1)/( a " T) correspond to ° aS = ° a • A comparison of these two sets of 
curves again verifies the statement following equation (20). 

In the following performance analysis of the smoother an important 
parameter is the optimum filter time constant ^p. This can be easily 
evaluated from the filter transfer function evaluated in Appendix C and is 
given by 


2T 


X F = ' 


Und-K^ 


( 21 ) 
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where K-| denotes the first component of the normalized Kalman gain BAKo* 
Figure 7 plots the normalized time constant (t f /t) as a function of 

2 2 

o 2 /a- 2 . For T = 0.1s. ( T p/T) varies between 2.5 and 14 as (o /a- ) 

a v r a v 

2 

varies over a large range. As is intuitively clear, higher values of cs & 

result in faster settling of the filter and this as shown later requires a 
smaller number of smoother delays to achieve the maximum possible smoother 

improvement. A comparison with figure 3 shows that ( T p/T) is of the order 
1/(2B p ). 

2. Optimal Smoother Performance with Linear Phase Detector 

In Figure 8 is plotted the smoother performance evaluated from 

2 2 

simulations as a function of the smoother delay and the ratio (a /a- ) 

3 V 

used in the smoother design, assuming a linear phase detector. The dotted 
curves in the figure plot the two-sided normalized loop noise bandwidth as 
computed from (16). As may be inferred from the figure, the two measures 
of performance are equal within the limits of statistical errors resulting 

due to finite data size. The minimum phase error variance (corresponding 

2 2 

to L = °°) varies over a range of about 0.3 to 1.4 for fa h- ) between 1 

a v 

to 100. It is also apparent from the figure that the amount of delay 
required to achieve asymptotic smoother performance has an inverse relation 
to (a a 2 /o- 2 ). In Figure 9 is plotted the real loop noise bandwidth 2B Lp>s 
as a function of normalized smoother delay (LT/Tp) f where is the time 
constant of the optimal filter. As is clear from the figure, the 

normalized value of delay required to achieve asymptotic smoother 

2 2 

performance does not depend significantly upon (a /a- ). 

3 V 

Figures 10 and 11 plot the results similar to those of figures 8 and 9 
respectively, when the actual process noise variance equals its design 
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2 


value, i.e., cr = c .A comparison of these two sets of figures shows 

a , S 3 

that the phase error variance can be at most 1.35 times more than for the 

2 


case of a „ = 0. For intermediate values of the noise variance, 

3 9 ^ 

2 2 

0 < a _ < a , the ratio would be smaller. 

3)0 a 

In figure 12, a comparison of smoother phase error variance is made 
for three different sampling periods T equal to .01, 0.1 and Is 
respectively. As is evident from the figure, whereas the optimum filter 
performance is dependent upon the sampling period, the asymptotic smoother 
performance depends only marginally on T. Moreover as the normalized 
smoother delay (to achieve asymptotic performance) is almost independent of 
T, the actual delay L has an inverse relation to T (from Figure 7 (t^/T) is 

r 

inversely related to T.) 


3. Smoother Performance with Nonlinear Phase Detector 

Figures 13-17 plot the performance of the smoother when the sinusoidal 
nonlinearity of the phase detector is included in the simulations. For 
these simulations, we consider a specific case wherein the design value of 

the process noise variance is equal to the measurement noise variance 

2 

a- (corresponding to a two-sided bandwidth of 1 Hz). From figure 13 it is 

2 

observed that when the simulated noise variance a c = 0, the effect of 

a,b 

nonlinearity is to degrade the normalized phase error variance Pg(1>1)/ 

2 2 

o- by at most a factor of 1.32 for cr- £ 2.2 (corresponding to the 

filter rms phase error of 27° for a linear detector and 30° for a nonlinear 
phase detector). 

In figure 14 is plotted the smoother performance with a nonlinear 

2 2 

phase detector for the case of o „ = o • Here the variation of 

a,S a 
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normalized phase estimation error variance depends much more strongly on 
o- . For a- < 1.1| (corresponding to a rms phase error of 28.5 at the 

phase detector output), the degradation is within a factor of 1.7 (2.3 dB). 

2 

The degradation can be much higher for larger values of o- . Figure 15 
plots the phase error variance normalized by the sampling period T vs the 
normalized smoother delay. Such a normalized result is expected to be 
approximately valid for different values of T. 

Figures 16 and 17 compare the filter and the limiting smoother 
performance as a function of optimal filter performance with a linear phase 
detector. Thus the inverse of the phase error variance of the filter with 

nonlinearity vs 10 log (P c /NqB l ) = 10 log (1 /c^ 2 ) is plotted in figure 16 

2 2 

for the case of o . = 0. Here, o, denotes the phase error variance of a 

a y & (p 

phase-locked loop with a linear phase detector and having a loop noise 
bandwidth of B L Hz and (P c /Nq) is the receiver carrier power to noise 

spectral density ratio. Figure 16 also plots the corresponding results for 
the smoother with and without nonlinearity. As may be inferred from the 

figure, for the case of a linear phase detector, the smoother provides an 
improvement of 5.6 dB over the filter. When the phase detector 
nonlinearity is taken into consideration, then for a range of 10 log 

(1 /n ) 2. 6 dB, the smoother still provides an improvement of at least 5.1 

<t> 

dB over the filter. Note, however, that the filter performance can itself 
be degraded by as much as 1.5 dB due to phase detector nonlinearity. Since 
the plots in figure 16 correspond to a fixed value of B L , these results 
indicate that with a smoother, the receiver can be operated with at least 
3.5 dB smaller (P c /Nq) ratio when it is desired to have a 0.1 or smaller 
value for the phase error variance. 


19 



Figure 17 plots the corresponding results for the case ° a g = °a * 
Here it is observed that the effects of nonlinearity are more dominant 
resulting in a threshold behavior in the smoother phase error variance. 

However, for 10 log (P C /N 0 B L ) 2 7.5 dB, the results in terms of smoother 

2 

performance are close to those for the case of & a g = 0. 

B. Frequency Tracking Performance 

In the following we present the performance of the filter and smoother 
in terms of the frequency estimation error. 

1. Optimal Filter Performance 

The frequency estimation error variance is presented in terms of the 
frequency tracking loop's normalized two-sided loop-noise bandwidth denoted 


by 2 Bf. This parameter can also be computed from (16) with G F>p replaced 
by Gp^j. - the second component of the transfer function matrix Gp(z). The 

actual two-sided loop bandwidth 28^^ equals (2Bf/T) or is equivalently 

2 

given by P p (2, 2)/(o- T). This parameter is plotted vs the design 

2 2 2 

parameter (o /a- ) for the case of o = 0 in figure 18. From these 

a v a , S 


graphics it 

is readily seen that approximately, 

2 \f 

= 0.32 ( a /a- ) 
a v y 

, T = Is , (a 2 /a- 2 ) < 

cl V 

2B Lf 

2 o-75 

a 0.35 (a /o-) 
a v 

; t = o.ls 

A comparison 

of (22) with (19,20) 

shows that approximately 

2B Lf 

SS 10 (2B ) 5 * 8 

Lp 

T = Is , (a 2 /a- 2 ) < 10 
a v 

2B Lf 

= .34 (2B ) 3,32 , 

Lp 

T = O.ls 

2B Lf 

S -30 (2B Lp > 3 , 

T = .01s 


( 22 ) 


( 23 ) 
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It is hardly surprising that for low values of T, the frequency estimation 
error variance is proportional to the third power of 2B p . Still more 
important is the fact that the frequency estimation error degrades at a 
much faster rate with B Lp as T is increased. 


2. Optimal Smoother Performance with Linear Phase Detector 

In figure 18 is also plotted the equivalent loop noise bandwidth of 

2 

the smoother 2B Lf s = Pg(2,2)/(a- T) for the case of Ts 0.1s. This result 
can be approximated by the following expression. 


2B 


Lf ,S 


.097 (a 2 /a- 2 ) 
a v 


0.72 


(24) 


A comparison with (22) shows that the ratio of B Lf to B Lf>s is given by 

3.6(a 2 /o- 2 ) 0 ’ 03 which for low values of ( a 2 /rj- 2 ) is approximately 5.6 dB. 

a v a v 

Figures 19 and 20 plot the normalized frequency estimation error 

variance for the smoother as a function of smoother delay corresponding to 

2 2 2 

o „ =0 and a = a » respectively. Comparison with figures 9 and 

a , b a , S a 

11 respectively shows that the smoother delay required to achieve the 
asymptotic improvement in frequency estimation is about two times that 
required for phase estimation alone. Even so, a normalized delay of less 
than two is adequate for realizing the asymptotic improvement. Also, as 
already discussed, the frequency estimation error variance is a much 

stronger function of (a /a- ) and the actual process noise variance a 2 . 

a v a, S 

3. Smoother Performance with Nonlinear Phase Detector 

In figure 21 is plotted the normalized frequency estimation error 
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2 2 

variance for the case of a = a- and the actual process noise variance 

a v 

2 2 

0 = 0. For o- 2.2, the filter and smoother performances are degraded 

by a factor of 1.18 and 1.24 respectively. Figure 22 plots the smoother 

2 2 

performance for the case of a _ = a 

a , S a 


C. Smoother Parameters 

Figures 23 and 24 plot both components of the asymptotic smoother gain 
K^OO as k -*• 00 (see equation (4)) for various values of the delay variable 
i, corresponding to two different values of the design parameter ( / a- 2 ). 
As is apparent from the figures, the smoother gain is negligibly small for 
i > 20 in both of these cases. 
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VII. Conclusions and Suggestions for Further Work 

The report has presented the performance of the suboptioal filter and 
smoother for the phase and frequency estimation of a sinusoid under the 
presence of both the process noise and observation noise. The performance 
predicted on the basis of linear estimation theory is in close proximity 
with the corresponding results obtained with simulations, when the phase 
detector is assumed linear. Such simulation results are applicable when 
the phase detector nonlinearity is taken into consideration and the 
receiver is operating under high signal- to- noise ratio conditions. Under 
such conditions, the smoother improves both the phase and frequency 
estimation errors as compared to the filter by about 5.6 dB. 

However, as the signal- to-noise ratio is reduced, the corresponding 
improvement is also reduced. Als<* such a difference is more severe when 
the process noise is present than when only the observation noise is 
present. In the presence of both the nonlinearity and the process noise, 
the smoother performance exhibits a marked threshold behavior, in that the 

smoother performance degrades rapidly with decrease in 10 log (1/o^) below 

4 > 

7.5 dB. However, for 10 log (P c /NgB L ) 2 7.5 dB, the results in terms of 

2 

smoother performance are close to those for the case of o, 5 = 0. 

From the simulation examples it also appears that the performance of 
the filter/smoother does not change very significantly if the actual 
process noise variance is smaller than its design value. 

Although not studied in this report, it is expected that the 
degradation caused by ignoring the process noise variance completely or by 
underestimating it in the design would be relatively much more severe. 
Thus in situations where the process noise variance is not known, an upper 
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bound for it should be used in the design. This however, can lead to much 
higher estimation error than the minimum possible value as it follows from 
the simulations (the actual phase error variance would be close to the 
value obtained when the actual process noise variance is equal to its 
bound) . 

Therefore, for the future research we propose an adaptive 
implementation which either explicitly or implicitly estimates the actual 
process noise variance. 

Acknowledgement: The author acknowledges a partial support provided by the 
California State University, Long Beach in the form of Faculty Research 
Award. 
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Appendix A: Smoother Equations 

In the following, the filter and smoother equations are presented 
which result by expanding the nonlinear measurement function h(x,k) into a 
Taylor series around the predicted estimate of x(k) and retaining only the 
linear and quadratic terms. Much of the simplicity in the resulting 
equations is due to the fact that 2^ c t and higher order harmonic terms are 
ignored. These harmonic terms result from the quadratic terms present in 
the filter/smoother equations. 

Denoting by h(x,k) the function A Sin( u) c t k + 3&"x(k)) and by h x (x,k) 
its partial derivative with respect to x, the smoother equations can be 
decomposed into two sets of equations termed filter equations and 
"smoother" equations as follows. Letting x^k) = x(k-i), x^k/j) denote 
the estimate of x^(k) on the basis of observations up to time j, one 
obtains [7-9]: 

Filter State Equations: 

x 0 (k+1/k+1) = $XQ(k/k) + Mg(k+1 )v(k+1 ) (A1) 

v(k+1 ) = y(k+1 ) - h(x(k+1/k), k+1) (A2) 

where M^k+I) is the filter gain. 

Smoother State Equations: 

x^k+l/k+l) = Xi_i (k/k) + M i (k+1)v(k+1), 0 < i iL (A3) 

where L denotes the smoother delay. 

Filter/ Smoother Gain: 

M i (k+1) = P i0 (k+1/k) h x o“ 1 (k+1), 0 i i i L (A4) 

Error Covariance Update: 

In (A4) P iQ (k/j) denotes E £{x(k-i) - x(k-i/j)} |x(k) - x(k/j) T }J 
with j = k-1. A recursive update for P i0 (k/k) is given as 
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P 10 (k+1/k + 1) = P 10 (k+1/k) - [P 10 (k+1/k)h^]a“ 1 (k+1) [P 00 (k+1/k)h;] T 

P i0 (k+1/k) = P i-1>0 (k/k)<f T , 0 < i £ L 

n n 

a(k) = h x P 00 (k/k - 1)h x + R(k) + T E E 

i,j=l r,s=l 

i,r s,j 

r 00 00 Sx* d-x? 3x 1 3xJ 


(A5) 

(A6) 


(A7) 


In equation (A7) x* denotes the ith component of x and denotes (i,j)th 
component of the matrix Pqq. The smoothing error covariance matrix P Ai 
(k/k) = E { x(k-i/k)x T (k-i/k) } , x(k-i/k) = x(k-i) - x(k-i/k), is given by 


P ±1 (k+l/k+l) - P i± (k+l/k) - P ±0 (k+l/k)h^ a _1 (k+l)h x P ±() (k+l/k);0 < i <L 

(A8) 


P ±i (fcfl/k) = P w w ( k/k) . 1 < i < L 

To specialize the algorithm to the signal model (1,2), h x and h xx are 
evaluated as, 


h -a/TaCos (0(k)); - A»6"M'Sin (0(k))*, 

x 3x 2 

0(k) = w t. + &A x(k/k-l) 
c k 


(A9) 


a _1 (k) = ( 2 P A 2 Cos 2 (G(k)) + R(k) + p 2 A 2 S in 2 (0(k)) 1 , P, = a'p oo * 

1 J (A10) 


The above expression for o(i)may easily be reorganized as 

o _1 (k)=A" 2 [ P(J) + hv l + R(k)]" 1 [l + d(k)Cos(2 0(k))l' 1 

(All) 

d(k) = ^ ; R(k) = R(k)/A 2 

P <P + % + R(k) 
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Expansion of the second term on the right hand side into a cosine Fourier 
series results in 


a 1 (k) * A 2 + h + R(k)J 1 |a Q + a^Cos(20(k))+. . . j 


(A12) 


where 

a 0 * 

a-|, etc., are the Fourier 

series 

coefficients. Thus 



J_ 1 

r T ' dx 


1 


a o 


TT J 

0 l+d(k)Cosx 


V 1 -d 2 

(k) 



2 

TT 

r Cosxdx 


2 

Jl - 1 , ..,1 

a l 


TT J 

0 l+d(k)Cosx 


d(k) 

{ V l-d Z (k) J 


(A13) 


Substitution of o"^(k) from (A12,A13) into (A4) results in 


M^k+1) = P i0 (k+l/k)AV2~ f [p^ + *5 P$ + R(k +1 )] 1 

A 2 | C^osCBCk)) + C^Cos (30(k) )+. . . j 


with 



(A14) 


Recognizing that P i0 (k+1/k) is baseband, the baseband part of the 
correction term in (A2,A3) may be written as 


MjCk+imk+1) = K^k+DnU+l) 

^(k+1) = A P i0 (k+1/k) C 1 [a 2 p + is a 2 p 2 + R(k+1)] 

/— 'r ^ $ 

n(k+1) = T2 y(k+1)Cos(0(k+1)) (A15) 

In (A15) K^k+1) is baseband, and the term M^k+1) Sin(©(k+1)) has been 
ignored as this does not contribute any baseband component. 

Similarly, the second term on the right hand side of (A8.) is given by 


3 i0 (k+l/k)U^P. 0 (k+l/k) A 2 [l + Cos (20(k+l) )J [ a 2 p^ + ± a 2 P(J) 2 + R( k+ 1 ) ] 

l 0 + a x Cos(20(k+1) ) +... } 
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and its baseband part is simply 

p 10 (k+1/k )( A 3^)( A 6WP i0 (k+l/k) |a 2 + h a 2 + R(k+1)] (a Q + a 1 /2) (A16) 
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Appendix B: Steady State Performance Equations for a Second Order Linear 

Fil ter/Smoother 

Denoting by R, the variance of the additive noise term — — v . (k+1) in 

y-j" i 

equation (8), the steady-state error covariance matrix Pp is the solution 
of the following algebraic Riccati equation. 


P = *Tp - p H(H'p h + R) h"p 1 §'+ Q 
P L p p p pJ w 

or 

P^ - P H(H"P H + R)”Vp = 4> _1 (P -Q)(<J>0 _1 
P P P p p v ' 


(B1 ) 


(B2) 


Denoting by P^j the ij th component of P p , the left hand side of (B2) 
denoted B equals, 


B = 


1 

(P n + R) 


P 11 R 


P 12 R 


P 12 R (P 11 P 22 



) + P 


22 



(B3) 


Letting A denote the matrix on the right hand side of (B2), one obtains 


A 11 = (P 11 - 2TP 12 + T V - Qn 5 
A 12 = (P 12 " TP 22 ) " Q 12 i Q i2 = Q 12 - t Q 22 
A 22 = P 22 " ^22 ; ^22 = Q 22 


^11 = Q 11 - 2t Q 12 + T Q 22 


Substitution of (B3) and (B4) results in the following set of equations in 
P 1 1 , P.J2 snd ^*22* 


(BA) 
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P 11 5 ■ (P 11 + *> | P 11 - 2TP 12 + t * P 22 - «ll) 

p 12 R- CP u +5) { p X 2 ~ T P 22 - ^12 } 

(P 11 P 22 " P 12 ) + P 22 P = (P 11 + P) <P 22 ~ 


(B5) 

(B6) 

(B7) 


Equation (B7) can also be written as 


or 


Q 22 P 11 " P 12 + R Q 22 ° 
P H + R - p^ 2 /Q 22 


(B8) 


(B 8) 


Substitution of (P^ + R) from (B"8) into (B6) and a little simplification 
yields 


P 


12 


T P 


22 




(B9) 


Substitution of (B9) and (B"8) in equation (B5) results in the following 
algebraic equation in P-|2* 


12 


<522 


Pl 2 - 2R + IQ i2 + Qi, 


12 


-2 

- T R P 12 + R =0 


(BIO) 


<22 


*22 


Defining P 12 = p i 2 / ^ Q 22 R * then substitution of ^11* °22» Q 12 in 
terms of Q 1 1 etc. from (BM) results in the following equation for P 12 , 


P't„ “ a p i2 ” (2 + b) P 12 “ 3 P 12 + 1 = ° 


12 




q 22 /r 


b = 


Q 11 “ TQ 12 


(B11) 
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In terms of the solution P^ 2 of (B11), the elements of the matrix P p 
can be computed as follows. 


12 = 

^22 R 


— r~2 

11 = 

R (P 12 

. 

^ Q 22 R 

22 

T 


12 


(B12) 




12 

Due to symmetry in the coefficients of the polynomial in equation 
(B11), it can be factorized as 

<? 12 - C 1 f 12 + l) <*12 - C 2 ? 12 + 

comparing the coefficients yields the solution for C 1 and C 2 as follows, 

2 ~ % + 4 (A + b)J 


( B1 3 ) 


Thus for the existence of the solution, one must satisfy the condition 

2 


+4 (4 + b) > 0 


Now P-j 2 is the solution of 


~ ? n ' c ? i2 + 1 ' 0 


(B14) 


(B15) 


Since it is required that |? 12 l > 1 (note that the solutions of (B15) are 
A and A" - ' for some real A) ; 

+\[c 


p 12 - * 


/ c 2 - A 


= h 


c - \jc 2 - 


c > o 


(B16) 


C < 0 


For the above solution to exist > 4. If only one solution of (B13) 
satisfies this condiition then the solution is uniquely determined. This 
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would be the case, for example, when 


|4 + b| < 4 or -8 < b < 0 (B17) 

If on the other hand both solutions of (B13) satisfy C 2 > 4, then the one 

satisfying P 2 2 > Q 22 ( p ^22 + ®22 = p? 22^ where superscripts F and P stand 
for filter and predictor respectively is the desired solution. 

With the knowledge of Pp, the steady-state filter error covariance 
matrix Pp is simply given as 

P p = p p - P p H (H'P p H + R)“lW p 

The smoother covariance matrix P iQ and the filter smoother gains can 
then be obtained from equations (C10-C13) of the Appendix C, which are time 
invariant versions of the relevant smoother equations. Moreover, the 
smoother performance for various values of the delay can be evaluated 
recursively from equation (7). Equivalently the asymptotic performance of 
the smoother can be evaluated from the following nonrecursive equation. 




( p “ P Q ) - P t'P 1 (P ~ PJ P' 1 ? P = P - P 
P S P P p s' p * p p r F 


where 


$ K q H' ; P s = lim lim P (k/k) J 

J L - * 00 k-» 00 ^ 


(B18) 


Letting, 


. 1 

P <TP = 
P P 


a Y 

6 6 
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<11 


then the solution for Pg = (P p - Pg) is given in terms of 


-p (P p - Pp) as follows, 


-1 


V 1 ’ 1 ) 


(l -ot ) -2oty -y 2 


V 1 .*) 

? s (1 ’ 2 > 

= 

-a8 l-(a5+8Y) -y6 


P F 0>2) 

? s (2,2) 


-S 2 -266 1-6 2 

” 


Pp(2,2) 





“ - 


(B19) 
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Appendix C: Transfer functions of the steady-state filter /smoother. 

Third-Order Filter : 

From equations (2), (3)» (8) and under the assumption of a linear 
phase detector, one obtains 


x(k+l/k+l) = <t x(k/k) + K q ABZ'^{ x(k) - x(k/k)} + K 0 A6?"w(k) 

+ K Q v^(k+l) 


(Cl) 


In equation (Cl), v - (k+1) = ==. v (k+1) , and {w(k} is a vector white 

V 2 1 

Gaussian noise process as appearing in equation (2). Subtraction of (Cl) 
from the second of equations (2) yields the following equation after some 
simplification and taking the Z transform on both sides of the equation. 


(zi - (i - ia") <j>)x f (z) = [i - iur]w(z) - K Q zv i (z) ( C2 ) 

In equation (C2) above, ic = BAK Q = [K 1 K 2 and X F (z),W(z) and V i (z) 

represent the Z transforms of x(k/k), w(k) and v^(k) respectively. Letting 
the vector function £ (z) represent the right hand side of (C2), the 
equation (C2) may be written as a set of following three equations in the 
unknowns X^z), X 2 (z) and Xg(z) (after substitution of $ from (14)). 


[z - (1-K ± )] x[ - (l-K^T X F 

K 2 X F + [z - (1-K 2 T)] X F 

- ~f — f r 

K. X. + K,T X + z 
3 1 3 2 L 


5 1 (2) 

(C3, 1) 

T(l-K 2 T/2) X F = C 2 (z) 

(C3,2) 

(l-K 3 T 2 /2)]x F = C 3 (z) 

(C3,3) 


The set of equations (C3) can be solved to obtain the following transfer 

~F 

function matrix between V i (z) and X (z) (setting W(z)=0) , 
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~F 

A8X (z) 
zV i (z) 


KjCz-l) 2 + (K 2 T + Js K 3 T 2 )(z-1) + K 3 T 2 
(z-1) [k 2 (z- 1 ) + K 3 t] 

C*-i > 2 k 3 


(C4) 


D(z) - (z-1) 3 + (Kj + K 2 T + K 3 -y-) (z-1) 2 + (K 2 T + K 3 T 2 ) (z-1) (C5) 


From (14) and (C5) it follows that 


+ K 3 T 


2 


Agl/3> X F (z) 
zV ± (z) 

~ 1 
D(z) 


(Kj+i^T-H^ -|-)(z-l) 2 + (K 2 T+ -y- K 3 T 2 )(z-l) 



(C6) 


From equation (8), in the absence of process noise w(k), the transfer 
function between n(k+1) and v^k+1), (or between n(k+1) and y(k+1)) is 
given by 


n(z) 


V.(z) 


Agr$ x F ( Z ) 

zV.(z) 

(z-l) 3 /D(z) 


(C7) 


By exploiting the linear property of the filter (or smoother) and 
considering the case when w(k) = 0 (x(k) m 0, implying that x(k/k) = 
-x(k/k)), one readily observes that the filter transfer function between 
y(k) and x(k/k) (or between v^(k) and x(k/k)) is given by 


A8G (z) £ AB JLli) 
* Y (z) 


2 

KjCz-l) + (K T+% K„T 2 )( z -i) + kt 2 


D(z) 


(Z_1) [ K 2^ z_1 ^ + K 3 t] 

(z-i ) 2 iL 


(C8) 
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The first term of the transfer function matrix on the right hand side of 
(C8) is simply A times the transfer function between the phase estimate and 
the input signal y(k). Similarly the second term represents A times the 
transfer function between the frequency estimate and the input and so on. 
Using frequency domain techniques, these transfer functions are used in the 
developed software to evaluate the filter loop noise bandwidth. 

Transfer Func tion of Third-Order Smoother: 

It is easily seen from equation (3) that the smoothed estimate (with 
smoother lag equal to L) of x(k+1-L) may be written as, 

x(k+l-L/k+l) = x(k+l-L/k+l-L) + JL (k+l)n(k+l) + 1 <k)n(k) 

(C9) 

+ ...+ KjCk+Z-DrHk+Z-L) 

The various smoother gains and the covariance matrices are now replaced by 
their respective steady-state values as follows. From equations (4) and 
(5) one obtains 

P i 0 (k+l/k) = P._ li 0 (k/k-l)r(k) (C io) 

where 

*00 = $(k) - AB$(k)K Q (k)r (C11) 

A nonrecursive solution of (CIO) under steady state can be obtained by 

~ ~ 
replacing $(k) and P 00 (k/k-1) by their steady-state values denoted by <p and 

Pp respectively. Thus, 
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(C12) 


lim p (k+l/k) = p ; l£i£L 

>. oo ‘ 

lim K (k+1) = ASP [7"1 1 £3 _1 

k- oo 1 P 


( Cl 3 ) 


Denoting by X s the Z transform of the smoother estimation error x(k+1- 
L/k+1), then the substitution of ( C 1 3 ) and taking the Z transform of the 
resulting error equation yield, 

_ A g jj z > = ab = A6 - z" (L_1) - (a6) 2 s _1 p [?"] n(z - ) 

V 1 (z) V i (z) V i (z) P V i (z) 

x | I + (J') _1 z -1 +...+(0" (L " 1) z _(L ‘ 1) } 


Substitution of (C6) and (C7) into (Cl 4) yields X^(z)/V^(z) or 
-X^(z)/V^(z) (taking x(k) = 0). Equation (C14) in its present form is more 
convenient for computer evaluation of the transfer functions and has been 
used for such purposes in the program SIM. However by summing up the terms 
in the finite series in (Cl 4), the second term on the right hand side of 
(Cl 4) may be written as, 

A6 z" (L " 1) (z-l)V 1 P $"(I-?"z) _1 { I -t^'3 L z L }it/D(z) (C15) 

In particular, the transfer function between the smoothed phase estimate 
0 S (k) and v^k) (ory(k)) obtained by premultiplying ( C 1 4 ) by %' is given 
by 


9 S (z) 

V ± (z) 


J_ z -(L-i) 
A 


{ic i (z-l) 2 + (K 2 T+5 s K 3 T 2 )(z- 1) + K-jT 2 
+ (z-1) V {I-? L Z L } {I-I'z} - Vk q }/D(z) 


(C16) 


Transfer functions to the smoothed estimates of the frequency and frequency 


(CIA) 

8 
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rate can be similarly obtained. 


First and Second O rder Filter/Smoothers; 

The results derived for the third order filter/smoother can be 
directly specialized to obtain the corresponding results for lower order 
estimators. For example, substitution of Kg = 0 in equations (C4), (C5) 
and considering only the first two components of X c result in 


AgX F (z) = -1 

zVAz) E(z) 


K z + (TK 2 -K 1 ) 
K 2 (z-1) 


E(z) 


2 

z 


(2-K 1 -K 2 T)z + (1-K 1 ) 


(C17) 


Similarly substitution of K 2 = 0 (and = K) in (C17) and considering 
only the first component of X result in the transfer function of the first 
ordeii f il ter. This after cancellation of the common (z-1) term may be 
written as, 

A X F (z) = _ K 

zV i (z) z-(l-K) (Cl 8) 


In the same manner the transfer functions rt(z)/V^(z) can be obtained 
from (C7) by substitution of Kg = 0 or K 2 = Kg = 0 respectively for the 
second order and first order filter respectively. Furthermore, as (Cl 4) is 
applicable to the smoother of any order, a minor modification of ( C 1 5 ) , 
( C 1 6 ) also yields the smoother transfer functions for lower order 
estimators. 
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Appendix D: Flow Chart of Simulation Program for Smoothers 

The flow chart of the computer simulation program SIM developed on the 
VAX system is given in Figure 25. The program contains both the time 
domain and frequency domain analysis of the smoothers. 



LOOP FILTER 


1 



x(k+iA+D 


x(k/k+l) 


x(k-l/k+l) 


x(k+l -L/k+l) 


gure 1. 


Smoother Implementation 


2 













LOOP FILTER 


Figure 2. An Equivalent and Simpler Implementation of the Smoother 
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Bandwidth of the Second Order Filter/Smoother vs (cr Jcr- 
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Figure 4. Real Two Sided Bandwidth and the Normalized Phase Error Variance vs (cr /cr— ) 
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Figure 6. Comparison of Filter Performance With and Without Process Noise 





Figure 8. Smoother Estimation Error Variance vs Smoother Delay 
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Figure 13. Smoother Phase Estimation Error Variance 

2 

vs L (Nonlinear Detector With a „ = 0) 
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Figure 15. Normalized Phase Estimation Error vs 

Normalized Delay (Nonlinear Phase Detector) 
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X NONLINEAR FILTER 
A NONLINEAR SMOOTHER 
O LINEAR SMOOTHER 


1 | | | I i i I 1 

5 6 7 8 9 10 11 12 

10 log (P c /N 0 B l ) — 


Figure 16. 
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NORM. FREQ EST ERROR VAR. P c (2,2)/(cr£ T) 



Figure 19. Normalized Frequency Estimation Error Variance vs 
Normalized Delay (<X c = *->) 
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Figure 22. Normalized Frequency Estimation Error Variance With 

2 2 

Nonlinear Phase Detector (cr „ = <r ) 
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Figure 23. Smoother Gain K.; l<i$L (<r /< T- = 1) 
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SMOOTHER GAIN K. 



Figure 24. Smoother Gain K.; l<i<L (a 2 /a- = 50) 
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Figure 25. Flow Chart of Simulation Program for Smoothers 
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Figure 25. Flow Chart of Simulation Program for Smoothers (Contd) 
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Figure 25. Flow Chart of Simulation Program for Smoothers (Contd) 
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